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Abstract 

We consider a holographic model of QCD at finite temperature with nonzero 
chemical potentials conjugate to -R-charge densities. A critical surface of the 
confinement-deconfinement phase transition is shown for five-dimensional charged 
black hole solution given by Behrnd, Cvetic and Sabra. On a special section of 
the parameter space, we find an critical curve being similar to the one expected in 
QCD. We calculate meson spectra and decay constants in the confinement phase 
of this section to see their temperature and chemical potential dependences. We 
could assure generalized GellMann-Oakes-Renner relation and the reduction of 
pion velocity near the critical point. 
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1 Introduction 



Recently, from the gravity/gauge correspondence [T], the probe brane approach to the 
gauge theory with flavor quarks has been largely developed in terms of the system 
of Dp/ D p+4: branes 0E1- Inspired by these works, phenomenological, 5d holographic 
models have been proposed to explain more quantitative properties of light mesons 
in a simple setting [H El IE] • And a simple model proposed in |U |3] is extended 
to the finite-temperature version which could cover both the confinement and the 
deconfinement phase. 

On the other hand, the thermal gauge theories with chemical potentials conjugate 
to i?-charges densities have been studied in the framework of the gauge/gravity corre- 
spondence PUJ-^T] . While the charge in such theories is not the one of QCD, it would 
be an interesting problem to apply our previous model [9 to these theories and to see 
the flavor meson properties in a thermal medium with chemical potentials conjugate 
to these charge densities. 

Among the models with i?-charges, we consider here the background configuration 
which has been found by Behrnd, Cvetic and Sabra jTU| as a solution to the equa- 
tions of motion of the five- dimensional M = 2 gauged supergravity. This solution is 
corresponding to the dimensional reduction of the spinning three-brane, and it is char- 
acterized by three chemical potentials and a horizon. The analyses are given here on a 
special section of these parameter space to study the dynamical role of these chemical 
potentials in the various quantities related to quark and meson. 

In the next section, we set our holographic model of i?-charge induced chemical 
potential. And a phase diagram, which is very similar to the one expected in QCD, 
of this model is given. In the section 3, various light meson properties are given and 
we find many interesting results, especially near the critical curve of confinement and 
deconfinement. And the summary is given in the final section. 



2 Model setting 

We use the .R-charged black hole solution to equations of motion of the five- dimensional 
J\f = 2 gauged supergravity found by Behrnd, Cvetic and Sabra [TU] . The background 
for R 3 three space is written as, 

7i- 2 ^ ( dz 2 \ 

dsl = ^^\-f 2 {z)dt 2 + n{dx i ) 2 + n- 1 



z 2 V J w v ' P{z) i 



f 2 {z) — H— (— ) 4 11(1 + «i) , K = H!H 2 H 3 , Hi = l + k^—) 2 , m = q t z 2 + , (2) 
z + t=i z + 

where (z = 1 ~ 3) are the charges of three Abelian gauge groups and z + represents 

the smallest value of f(z) = 0. The radius of AdS 5 is taken as unit. The Hawking 
temperature of the background is given by 

2 + Ki + K 2 + K-3, ~ KlK 2 ^3 rr, / Q \ 

1 — ; -tO W) 

2J(l + « 1 )(l + « 2 )(l + « 3 ) 



1 



where To = 1/(ttz+) denotes the Hawking temperature in the case of the zero chemical 
potential. Here the chemical potential is defined through the three gauge potentials, 
A\(z), of the gauged supergravity as 



Hi = A\{z) 



Z=ZJ r 



Z + (l + Ki 



n( i +^ 



,1/2 



(4) 



1=1 



Since the bulk configuration depends on the gauge potential, we should consider 
the thermal properties of the system by considering the Gibbs potential in which the 
chemical potential is fixed. The density of the Gibbs potential Vq is given as 

We find from Vq the condition of thermodynamic stability by demanding that the 
system is on the local minimum point of G in T, fii space. This requirement implies 
the following constraint * on Ki 

2 - Ki - K 2 - K 3 + Ki K 2 K 3 > . (6) 

Although various analyses for this system have been given from the viewpoint of ther- 
modynamics, we consider only the stable region, which satisfies (0), in the parameter 
space. And we give a critical surface for quark confinement and deconfinement phase 
transition from a different point of view based on an approximate holographic model. 

Our idea is as follows. Usually, the 5d black hole solutions are used to study the high 
temperature phase of the dual 4d gauge theories without quark confinement. However, 
for the same black hole solutions, the low temperature phase with confinement can be 
realized by introducing an infrared cutoff for z, z m , such as z m < z + . In this case, 
the physical region is restricted to the region of < z < z m . As a result we can see 
discrete meson spectra, which are characterized by z m , and their properties at finite 
temperature. 

On the other hand, for the case of z m > z + , the physical region is restricted to 
< z < z + and we find the usual deconfinement background. Then, in this formalism, 
the critical temperature T c for the confinement- deconfinement phase transition is given 
as T c = l/(7rz m ) for the case without the chemical potential. The value of z m is 
determined by using some experimental data of meson mass, for example by the vector 
meson (p meson) mass, and we obtain T c ~ 100 MeV. 

When the chemical potentials are introduced by the model given above, we obtain 
in a similar way the critical surface as 

T _ 2 + Ki + K 2 + K 3 - K!K 2 K 3 1 

2J(1 + /C 1 )(1 + K 2 )(1 + /C 3 ) ^ ' 



This is considered in the 4d parameter space of temperature and three chemical po- 
tentials. In order to visualize it, we give a typical example of such a critical surface in 
a 3d parameter space for K\ = k 2 in the Fig^ 

*We notice that the sign of the last term in Eq. © is opposite to the one of This point is 
however important. 
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Fig. 1: (A) Critical surface in the space of ((«i = k<i —)k, k 3 ,T). The surface on 
the hill represents the thermodynamic stable region. (B) Critical line obtained for 
(Hi = ^2 = ^3=)/^ an d T. This is the diagonal section of the critical surface in the 
Fig. (A). T and \x are scaled by l/^z m and z m respectively. 



The surface on the hill in Fig.QjA) represents the critical surface on which the black 
hole solution is thermodynamically stable. Outside of this hill, the bottom region of 
the valley, the i?-charged black hole is considered to be unstable. So we must consider 
some other background in this region. We fortunately find an interesting slice given by 
the section of Ki = k,2 = ^3 = ft, from (T, n) = (1, 0) to (0, 2), on the critical surface. 
This section provides a phase diagram similar to the one expected in QCD. It is shown 
in the Fig.^B). 

Hereafter, we make analyses on the special section, and the k (/i) and the T depen- 
dence of various physical quantities are studied. 



3 Meson spectra 

We consider the 5d meson-action proposed in [H E] but use the background ((TJ) to treat 
the system of finite T and k. The meson action is 



'S'meson — / d xdz a/ ~Q Tr 



1 ;(L MN L MN + RmnR MN ) - \D M ^ - Mlm 



where the case of two flavors (Nf = 2) is considered; notations are the same as in j^j. 
This action is considered as the quark part represented by the probe brane in the string 
theory. And the fields on the brane is considered to represent meson states composed 
of a quark-antiquark pair. 



3.1 Scalar field 

The scalar is defined as $ = Se^ 1 '"' and v(z) = 2(S). The 5d mass of $ is M| = —3, 
since the bulk field corresponds, in the gauge theory side, to an operator with the 



conformal dimension A = 3. The equation of motion for v is obtained as 



Z 2 f 2 



v(z) = 0, 



(9) 



and v is numerically evaluated in this paper. The asymptotic form near z = is 
v(z) ~ m q z + cz 3 , since the background (jTJ) tends to the AdS background there. Two 
constants, m q and c, are identified with the quark mass and the chiral condensate, 
respectively. 

The v diverges logarithmically at z = z + because of the factor /. We have confirmed 
it explicitly for the case of T > and k = [5| with the analytic solution of v(z), and 
numerically for the case of finite k. However, we can evade the dangerous region z ~ z + , 
when z + is larger than z m to some extent. The present analysis is then focused on the 
case z + > 1.1 z m for which the logarithmic divergence is not realized in the physical 
region z < z < z m . In analyses shown below, we fix T and see the k dependence of 
physical quantities, since the T dependence has already been discussed in the case of 
k = [9 j. So the case z + > 1.1 z m is realized in the lower k region, say k < 0.9 n c (T), 
where the critical density k c (T) is a function of T and satisfies the condition z + = z m 
for each T, that is 



2 + 3k c (T) - K c {Tf 1 
2(1 + k c (T)) 3 /2 ttT 



(10) 



The treatment of the dangerous region near k = k c (T) will be discussed later. 

In the present framework, without any inconsistency, we can allow parameters m q 
and c to depend on T and k. However, we do not find a way of determining the T and 
the k dependence, although in real QCD we know that the chiral condensate c does 
depend on T. This point is the defect of the present model. In the present analysis, 
we simply assume m q and c are constant in the region k < 0.9 k c (T). This would be a 
reasonable assumption for m q , but it is just a simplification for c 

The fluctuation of S, which is defined as S = v(z)/2 + a, can be observed as a 
singlet meson state (a). Here and hereafter we consider the static mode, di<p = 0, for 
any field in order to derive the mass in a simple way. So the invariant mass, or pole 



mass, is defined here as — <9 t 2 ( 
4d mass term to Eq. as 



m 



Then the equation for a is given by adding this 



3 f2 



d z + 



m 1/3 h 



z 2 f 



H -m 

f 4 



a 



This equation is independent of v(z). The improvement of this defect is postponed to 
the future. 

The discrete mass-spectrum is obtained by solving this equation with the boundary 
conditions, a(z)\ Zo = d z a(z)\ Zm = 0, where z is the UV cutoff which is taken to zero 
after all. The mass depends on z m , T and k. First, we determine the value of z m so as 
to reproduce p meson mass at T = k = 0. The resultant value is l/z m = 0.323 GeV 
[7]. It is then possible to estimate the k dependence of a meson mass. 

The k dependence of m is predictable from Eq.(jTTJ). In the equation, the mass 
term m 2 Ti,/f A will have a strong k dependence through the factor f A , if m has no k 
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Fig. 2: The k dependence of f 2 /\^H with T fixed. Three curves correspond to three 
cases of T = 1, 40, 80 MeV, respectively. Note that k c = 1.989 for the case of 
T = 1 MeV, 1.509 for T = 40 MeV and 0.858 for T = 80 MeV. 



dependence. However, the mass term should have a weaker k dependence, since so 
do the other terms. Hence, the strong k dependence of f A is suppressed by that of m 
in the mass term, indicating that the k dependence of m becomes similar to that of 

f 2 /Vn. 

Fig. |21 shows the k dependence of the factor f 2 /\^H with T fixed. Three curves 
correspond to three cases of T = 1, 40, 80 MeV, respectively. The k varies in a region 
< k < k c . The phase transition takes place when k = k c . For lower T cases such as 
T = 40 and 1 MeV, as k increases, the factor increases once and finally tends to zero 
in the limit of k = k c . 

The k dependence of p/y/TL is compared with that of m calculated numerically. 
The latter is shown in Fig. |3] for the case of T = 40 MeV. The k dependences of m is 
found to be similar to that of f 2 /VTi. This similarity is true for other T. Thus, we 
can see a tendency that m decreases as k goes to k c . 

For a meson, it is possible to take the Dirichlet condition a(z)\ Zm = 0. It makes 
m larger by ~ 50%, but the ^-dependence of m is similar to the case of the Neumann 
condition. 



3.2 Vector meson 



The gauge bosons are separated to the vector and the axial vector, Vm and Am, and 
are defined as L M = Vm + A m and Rm = Vm — A M , respectively. 

First of all, we consider the vector mesons. The linearized equation for the spatial 
component Vi is given as 



m 2 H ^ z 



P 



-d z 



■d z 



V, = 



(12) 



gauge is taken. We adopt the same boundary condition Vi(z)| Zo = 
Zm = as the case of the scalar meson. This equation yields the discrete eigen- 



where V z 

d z Vi(z]\ 

values m 2 = w? n under the boundary condition. In order to see the decay constants, it 



is convenient to expand Vi(x, z) as Vi(x, z) = J2 n vl n \x)h v (z) for each mass eigenstate. 
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Fig. 3: The k dependence of masses and decay constants with T = 40 MeV fixed. In 
(A) the curves show m^, m p , m a and m ai from the bottom. In (B), the curves denote 
F^/2, F*, sfFp an d \J F ai from the bottom. The F\ increases as k goes to k c , while F* 
decreases. Other parameters are fixed as m q = 2.26 MeV and c = (0.333GeV) 3 . 



The spatial component of the wave-function for each mode is normalized as 



m dz^—{hl{z)f = \. (13) 

20 Z J 



The factor Ti. 2 ^ 3 /(zf 2 (z)) in the integrand is important in the sense that it determines 
the k dependence of the decay constant F p . 

Equation (fT2"|) is also independent of v(z), and the K-dependence of the lowest mass, 
which is identified with p meson mass, is shown in Fig. El The p-meson mass vanishes 
or becomes small as k goes to k c because of the mass term m 2 7i/ f A . 



3.3 Axial-vector meson and 7r-meson 



The linearized equations of motion for the axial vector and the pion 7r are obtained 

as 



. p 



+ d 2 + 



n 1/3 n/ p 



P 



-d 



•a 



9l 



o z ( ) • d 9 . 



H 2/3 



z 



9 5 - 



,H 1/3 v 2 



z 2 P 



z 2 P 
A ± = 0, 



At 



d 2 + 



OA • d 7 _ 



H 2/3 



z 



,H 1/3 v 2 



z 2 P 



(tt + ^) = 0, 



v 2 f 2 

2 q 2 u J 



d z % = 0, 



0. 



(14) 
(15) 
(16) 
(17) 



where is decomposed into the transverse and the longitudinal part, A^ = A^+d^ip, 
and A z = gauge is taken. For simplicity, flavor index is neglected and gc, is determined 
from the vector current two-point function at T = k = jl]. These equations are 
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solved numerically under the boundary conditions, A^±(zo) = d z A /J> ±(z m ) = and 
ip(z ) = d z tp(z m ) = tt(z ) = 0. 

Differently from the case of vector and a mesons, the quantities of axial-vector 
and pion depend on five parameters m q , c, z m , T and k through v(z), f(z) and H as 
seen in Eqs. (|14|) -(|17 |) . For the consistency between the vector and axial- vector meson 
sectors, here we take the same z m as in the vector meson sector. Parameters, m q and 
c, are determined to reproduce the experimental values, rh^ and F n , of m n and 
at T = k = 0; the resultant values are m q = 2.26 MeV and c = (0.333 GeV) 3 . This 
parameter set reproduces masses and decay constants of ir, a, p, a\ mesons at T = k = 
within ~ 10% error [H Ej. The present model is also successful in reproducing the 
T dependences of masses and decay constants of n and p mesons calculated with 
the mo del- independent methods such as lattice QCD calculation [20] and the chiral 
perturbation theory [T9] . 

In principle, the chiral condensate c depends on T and k, but we simply assume 
that c is independent of T and k. As shown in Fig. El masses of axial-vector meson 
(ai) and pion also become small as k goes to k c . 



3.4 Decay constants 

The decay constants are obtained from the wave functions as |3J 



gl 



\d 2 K 




_ dz 2 


zo- 



1 2 



1 d z A, 



(o) 
o.u 



gl 



where Ai±(x, z) = 
Furthermore: i»> 

|(0) 



En>l 
|(0)n 



at\x)K{z) 



and h^(z) is normalized as h%{z 



(o) 



) given in (JT3J). 

is the solution to Eq. (fTHjl (Eq. (fTlj) with = 0), satisfying 
Aq'1_ l (zo) = 1 and d z A^- ± (z m ) = 0, and F^ ,s are the timelike and spatial components of 
the pion decay constant. Obviously, F\ and F^ are different from each other at finite 
T or k, while = at T = k = 0. 

We can see from Figs. El that the pion decay constants decreases as k increases. 
Eventually, we find m(ft)/m(0) ~ F(k)/F(0). Note that F^ goes up slowly as k 
increases. This is an exception and important in determining the k dependence of pion 
velocity mentioned later. 



3.5 GOR relation and pion velocity 

In real QCD, a generalized GOR relation, m^F* 2 = 2m q c, is expected for finite T [^TJ 
IT^j . In the previous work |9,, we confirmed that the relation is realized in our model 
for the case of T > and k = 0. It is then of interest whether the relation persists even 
for finite k. This is tested by calculating m v and F* numerically. Fig. 0] shows the m q 
dependence of m 2 and 2m q c/ F^ 2 by the solid line and the closed circles, respectively. 
Here the case of T = 40 MeV and k = 0.8 k c is taken as an example. The m 2 (T) shown 
by the solid line tends to zero as m q decreases, as a reflection of the Nambu-Goldstone 
theorem at finite T and k. Comparing the solid line and the closed circles, one can see 
that the generalized GOR relation is satisfied. 

The pion velocity v n in the thermal midium with finite p can be estimated ana- 
lytically in the chiral limit, m q = 0, from F* and as v n = F°/F* j2U E2j- In the 
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Fig. 4: The m q dependence of m 2 . Other parameters are fixed at T = 40 MeV, 
k = 0.8 k c , k c = 1-509 and c = (0.333 GeV) 3 . The solid line corresponds to m 2 (K) and 
the closed circles to 2m q c/ F^{k). 



outer region of the critical line, i.e. when k > K c , we expect the restoration of chiral 
symmetry. So we can set as v(z) = for any k near k c , at least when the transition is 
the second order. This v has no divergence, even when z + = z m . We can then consider 
the limit k — > K C (T). The solutions of Eqs.(JHJ) and (fTSJl are given, under the boundary 
conditions v4j i0 | 2 =o = 1, as 

A, = l + h d Zl ^—, (18) 
Jo f 2 

A = l + b £dz 1 ^ JE . (19) 

The constants 6 ,i are determined by the boundary conditions, d z Ai$\ z=Zm = e, and the 
limit e — > is taken after obtaining v n . We then get 

TL{z m ) z + 

We are now considering the limit of z + — > z m , so v% approaches to zero linearly with 
respect to the difference z + — z m . This implies v^ 2 oc 1 — jr for the path of z + — > z m 
with k fixed and that v n 2 oc k c — k for the path with T fixed. Then we can assure the 
critical exponent v — 1 for both cases. 

Finally, in the chiral limit, we discuss the behavior of the chiral condensate c near 
k = k c (T). When c is finite, v(z) diverges at z = z + , although the divergence is loga- 
rithmic and then very weak. When k = k c (T), the divergence is realized in the present 
model which considers the region zq < z < z m . As mentioned above, the k depen- 
dence of c is not determined within the present model. However, if the deconfinement 
transition and the chiral one take place simultaneously at k — k c (T) and the chiral 
transition is the second order, we might evade the logarithmic divergence of v(z), since 
c is considered to tend to zero more strongly as k approaches K C (T). In this situation, 
we can evaluate v n at k = k c (T). 
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4 Summary 



In this paper, we constructed a holographic model of gauge theory at finite temperature 
with non-zero chemical potential (//) conjugate to i?-charge denisty (k), introducing 
the five-dimensional i?-charged black hole solution of Behrnd, Cvetic and Sabra. The 
gauge theory thus constructed has two phases, confinement and deconfinement ones. 
The phase diagram of the gauge theory is richer than that of QCD in the sense that 
the gauge theory has three chemical potentials conjugate to i?-charge densities while 
QCD has only one chemical potential conjugate to the baryon number density. In the 
three-dimensional phase diagram of the gauge theory, a special section is similar to the 
one expected in QCD. 

The present model makes it possible to evaluate meson masses and decay constants 
nonperturbatively even in the thermal system with non-zero chemical potential. We 
then calculated these quantities in the confinement phase of the special section and 
found that these quantities decrease as .R-charge density k increases with T fixed. 
Furthermore, we confirmed that the Nambu-Goldstone theorem and the generalized 
GellMann-Oakes-Renner relation persist in the thermal system with finite chemical 
potential. 

We estimated the pion velocity as v w = on the critical line of the phase diagram, 
assuming that the deconfinement phase transition and the second-order chiral one 
take place simultaneously. The present result is consistent with the result of Son and 
Stephanov [221 based on the sigma model, but somewhat deviates from the result of 
Ref. j2H|. However, the present result supports the statement that the measured pion 
velocity 0.65 [23], deduced from the pion spectra observed by STAR [IB] at RHIC, 
would be a signal of QCD phase transition. 

The present model has two defects. One is that the profile function v(z) diverges 
on the critical line, although the divergence is logarithmic and then very weak. The 
other is that the T and the k dependence of c are not determined within the present 
model. As for the chiral limit, however, if the deconfinement phase transition and the 
second-order chiral one take place simultaneously, we can expect that c tends to zero 
as T and k approach the critical line. For such a vanishing c, we can expect that v(z) 
does not diverge on the critical line and then that the present model becomes reliable 
even near the line. So it is quite interesting to find a way of determining the T and 
the k dependence of c in the present holographic model or its extension. 

Finally we comment on the bold procedure, the hard IR cutoff, which is essential 
in the present model. This cutoff seems to be bold, but we could understand it very 
naturally. The action S^son, Eq.(jSJ), is regarded as the effective 5d action for the probe 
brane(s). In the phase of chiral symmetry breaking, we find a minimum value of r, 
r m , for the probe brane by solving its profile function f|, so the theory on the brane 
should be restricted to the region r m < r. This r m could be identified with the IR 
cutoff introduced in the present model. 
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